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t-H Abstract 
t-H ■ 

C^) 'Low-energy structures of a hybrid system consisting of a polymer and an attractive nanowire substrate as well as the 
C^l thermodynamics of the adsorption transition are studied by means of Monte Carlo computer simulations. Depending on 
i— ( Structural and energetic properties of the substrate, we find different adsorbed polymer conformations, amongst which 
are spherical droplets attached to the wire and monolayer tubes surrounding it. We identify adsorption temperatures 
and the type of the transition between adsorbed and desorbed structures depending on the substrate attraction strength. 
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1. Introduction 

In a systematic study of low-energy states of polymers 
adsorbed to a stringlike substrate, we recently found a vari- 
ety of different conformational phases by just changing two 
basic substrate properties [l[ . Among these structures arc 
spherical droplets attached to the string as well as barrel- 
like monolayer conformations surrounding the string. The 
latter conformations exhibit similarities to single-walled 
carbon nanotubes. 

Generally, the adsorption of polymers on material sur- 
faces or substrates is a crucial and nontrivial process in 
nature and nanotcchnology. It is also known, for exam- 
ple, that the adsorption process or, in fact, the poten- 
tial of a polymer to adsorb at a semiconductor surface 
depends essentially on details like the exact position of 
single monomers in the primary structure of a hetero- 
polymer Q. However, specific hybrid systems composed 
of inorganic matter and polymers potentially facilitate the 
development of completely new nanotcchnological devices 
like sensors for specific single molecules or devices for ul- 
trafast photonics 0, 0] ■ 

Hence, fundamental investigations are inevitable for 
'a better understanding of such systems. Computational 
studies of coarse-grained systems have proven to be ad- 
equate and quite useful for this purpose in recent years, 
both for predicting and interpreting specific and basic be- 
havior of polymer adsorption 0, ■ 

Particularly interesting are systems with a cylindri- 
cal substrate, like carbon nanotubes. Many of the spe- 
cial properties of these structures, that make them poten- 
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tial candidates for technological applications, can be con- 
trolled, influenced or amplified by coating them with poly- 
meric material 0,0]. In previous computational works, for 
example, the wetting of cylindrical fibers with polymers or 
the helical wrapping of single polymers around nanocylin- 
ders were studied j^, E3] ■ 

In a recent study, we have revealed a more general 
picture of the adsorption of polymers at different ultra- 
thin cylindrical substrates jl[ by performing generalized- 
ensemble Monte Carlo simulations [lll - [l3l | for a general 
coarse-grained model of this hybrid system. Here, after 
describing technical details (Sects. [2] and [3]), we focus on 
specific conformational transitions at very low tempera- 
tures (Sect. [3}, and discuss in a thermodynamic analysis 
the adsorption transition at comparatively high tempera- 
tures (Sect. [5]). 

2. The model 

For our adsorption study, we employ a coarse-grained 
off-lattice model, where the polymer consists of identical 
monomers which are represented by beads without internal 
structure. These arc connected sequentially by stiff bonds 
of unity length. In order to facilitate future enhancements 
and the comparison with previous studies (see, e.g., Q), 
we introduce a weak stiffness between the bonds, i.e., the 
polymer is not strictly flexible. The polymer is placed into 
a simulation box which also contains an attractive thin 
string located in its center. Its orientation defines the z- 
axis. The edge lengths of this box in x and y directions 
with periodic boundary conditions are chosen to be twice 
as large as the length of the completely stretched polymer. 
We note that the polymer is not grafted to the string and 
may move freely in space. 
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The total energy E of the polymer consists of contribu- 
tions from the Lennard- Jones interaction Vlj between all 
pairs of nonadjacent monomers, a weak bending stiffness 
Vbend and the monomer-string interaction V s tnng- 
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where nj is the geometrical distance between two mono- 
mers i and j, Qi is the angle between the two bonds con- 
nected to monomer i, and r z -i is the distance of the ith 
monomer perpendicular to the string. The parameters are 
set as follows: e m = er m = 1, such that Vlj(2 1 ' 6 ) = — 1. 
The bendin g st iffness is chosen to be comparatively weak, 
k = 0.25 0,113. 

The interaction V^tring between the string and the mono- 
mers is also based on a simple Lennard- Jones potential, 
where the wire is assumed to have a homogeneous "charge" 
distribution The string potential can then be con- 

sidered as the limiting case of the potential of a cylinder Q 
in the limit of vanishing radius and keeping the overall 
charge fixed [lij . 

Alternatively, the Lennard- Jones potential for the in- 
teraction between a monomer and the string can be inte- 
grated out along the string axis to yield (|Tc|) . 
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where <7f is the van der Waals radius of the string and can 
be considered as its effective "thickness". It is related to 
the minimum distance r™ m of the string potential via 
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For convenience, we scale the string "charge" density r]{ 
in such a way that the minimum value of the potential is 
String (»"™ in ) = —1 independently of at, i.e., we set ijf « 
0.53 aj 1 . Fig. Q] shows the correspondingly scaled string 
potential for different values of (Tf. 



3. Simulational details 

Polymer systems are known to possess highly nontriv- 
ial, rugged free-energy landscapes [l6| . For the simulation, 
we have, therefore, applied generalized-ensemble Monte 
Carlo methods. The ground-state energies have been esti- 
mated by using various, but conceptually similar, stochas- 
tic methods such as parallel tempering, Wang-Landau, 



and multicanonical sampling [ 1 lM 1 3t Il7j . as well as es- 
pecially designed optimizing approaches like energy land- 
scape paving, where the energy landscape is deformed ir- 
reversibly during the simulation [HS|. The efficiency of 
all stochastic methods strongly depends on the conforma- 
tional update set used (and, of course, on the fine-tuning 
of each method). If a move set is chosen reasonably well, 
all methods lead to comparable results in similar times. 
In addition, we have refined low-energy states by standard 
deterministic optimization techniques such as the conju- 
gate gradient method. 

For the estimation of the density of states and hence 
all thermodynamic quantities, we employed the Wang- 
Landau method lH 21 1 for the determination of the mul- 
ticanonical weights and performed a final production run 



in the multicanonical ensemble [ill, [l2j, |17[ . Independently 
of the values of the potential parameters, we partition the 
simulated energy interval into 10 000 bins in each simu- 
lation. The actual bin size hence depends on the energy 
range delimited by the putative ground-state energy and 
a fixed boundary on high energies. 

As known and mentioned above, the choice of the up- 
date scheme is crucial for the efficiency of any Monte Carlo 
simulation in general. For conformational changes, we ap- 
ply here a variety of update moves, see Fig. [2j includ- 
ing local crankshaft [Fig. Efa)] and slithering-snake moves 
[Fig. |2fb)], as well as global spherical-cap [Fig. [He)] and 
translation moves. Sets including these steps have been 
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Figure 1: The interaction potential between a monomer and the 
string. Note that by scaling the string "charge" density with Uj~ , 
the minimum value of the potential is —1 independently of Of. 



Figure 2: Conformational update moves used in our simulations, 
(a) local crankshaft move, (b) slithering snake or reptation move, 
and (c) global spherical cap update. 



2 



found to work quite well in previous studies as well [22H24j . 
The crankshaft move is just a rotation of a single monomer 
around the virtual bond between its neighbors, or, if the 
monomer is an end-monomer, around the elongation of 
the one neighboring bond. For the slithering-snake up- 
date, we cut a monomer from one end and paste it at the 
other end keeping the bond vector fixed. Both updates 
induce small conformational changes of the whole chain, 
whereas the latter one enables the polymer to leave very 
dense, adsorbed conformations. 

The spherical-cap update consists of the shift of 1 < 
n < N monomers by a small constant vector keeping the 
bond lengths at the nth monomer fixed. It hence allows 
for larger steps in the conformational space compared to 
the former local updates. The global translation update 
finally allows for a direct displacement of the chain relative 
to the string, which in the entire simulation remains fixed 
in the box. In any Monte Carlo step, we choose the dif- 
ferent moves randomly with equal weight. We convinced 
ourselves that this provides a reasonable sampling of the 
conformational space of this model system on large scales 
as well as locally. It behaves in general not worse and 
for the ground-state search even better than a procedure 
consisting of a few global moves followed by much more 
local moves, whereas the ratio depends on the system size. 
However, such a procedure has been found to be favorable 
for other problems [24[. 

Alternative, more sophisticated update moves like bond- 
bridging or monomer-jump moves |23l426l | have not been 
included in our update set as they are rather time con- 
suming and would apparently not improve the principal 
findings of the present study. Though, they are necessary 
for studying structural transitions in the dense and crys- 
talline regime 22, [23|, or for the investigation of much 
larger systems. In the following, we discuss the structural 
properties of a polymer with N = 100 monomers. 

4. Ground states of the system 

There are two parameters in the interaction potential 
between the string and monomers [Eq. (Tic)) and Fig.Q], the 
van der Waals radius proportional to the effective thickness 
of the string, Of , and the adsorption strength ef. By vary- 
ing these two parameters, we have recently constructed 
the complete conformational phase diagram of low-energy 
structures [l], [2?]]. By introducing suitable observables, 
four major structural phases have been identified. For 
small values of ef and at, i.e., for very thin and weakly 
attracting strings, we find globular or spherical polymer 
droplets inclosing the string (phase Gi). In this case, the 
polymer structures are similar to that in bulk under poor 
solvent conditions. The string docs not influence the shape 
of these structures but affects only the internal structure 
of the droplet. 

Fi gurc [3] shows conformations with ef — - e m — 1. In 
Figs. [3Ja) and[3£b), droplets inclosing the string are vi- 
sualized. When increasing the van der Waals radius of 




Figure 3: Low-energy conformations for ef = 1 and (a) erf = 0.55, 
(b) 0.573, (c) 0.647, (d) 1.0, and (c) 1.5. (a) and (b) correspond to 
phase Gi, conformations in (c)-(e) belong to phase Ge. 



the string, monomer-monomer bonds inside the droplet 
will be broken and, hence, the string is excluded from the 
droplet but still it is attached to it (phase Ge). The ra- 
dius at which this rearrangement occurs depends on ef. 
For ef < e m , i.e., where the string attraction is not signif- 
icantly stronger than the monomer-monomer attraction, 
the transition occurs, roughly, when the diameter of the 
string becomes comparable to the equilibrium distance be- 
tween two monomers. When further increasing the string 
radius, the string moves outward and the structure ap- 
proaches the bulk conformation. See Figs. [3tc)-(e) for 
examples of conformations at different <jf values, thus vi- 
sualizing the described "process" . 

Increasing at any given value of <7f the string attraction 
strength ef, globular structures (Gi, Ge) start to deform 
and to lose spherical symmetry at ef > 3. For erf < 1.5, we 
observe a transition from phase Gi directly to the barrel 
phase B, which is characterized by closed, stretched con- 
formations with cylinderlike shape wrapping around the 
string. In the extreme case of very strong string attraction, 
polymers form monolayer tubes. For Of > 1.5, the poly- 
mer first adopts clamshell-like conformations (phase C) 
before it forms barrel structures in phase B. Interestingly, 
the evolution from spherical droplets to monolayer tubes 
involves the formation of distinguishable monomer layers. 

For illustration, we plot in the upper row in Fig. U 
the radial distribution of monomers with respect to the 
string of certain low-energy structures. One finds accumu- 
lations of monomers at different distances from the string, 
i.e., in different layers. The position of the first layer cor- 
responds to the van der Waals radius of the string [cp. 
Eq. ([3]), <jf = 5/3], whereas the location of the higher- 
order layers is connected to the equilibrium distance be- 
tween the monomers, corresponding to a m . In the lower 
row, respective structures are depicted. In Fig. Eta), we 
plot the radial distribution of monomers in a conformation 
from phase Ge with ef = 1.5. The emergence of different 
peaks in that function can be observed. A clear 3-layer 
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Figure 4: Top: Radial distribution functions of low-energy states with <Tf = 5/3 and (a) ef = 1.5, (b) 3.0, (c) 3.5, and (d) 4.5. Bottom: 
Visualizations of the respective structures. Different colors or shapes encode different monomer layers, i.e., regions within certain distances 
to the string. 



structure can be identified in Fig. Etb), where a typical 
conformation in phase C is shown (ef = 3). In Fig. SJc), 
a two-layer barrel-shaped conformation is shown (ef = 3.5) 
which transforms into a monolayer tube at ef > 4.5, as de- 
picted in Fig. HJd). We would like to note here that, in 
particular, these monolayer tubes exhibit interesting sim- 
ilarities to other structures in nature like, for example, 
carbon nanotubes 111 r27H29f . 



5. Thermodynamics of the adsorption 

Finally, we short comment on thermodynamic proper- 
ties of the adsorption transition [l5|. For estimating the 
finite-system transition temperature, we first identify peak 
positions in the canonical specific-heat curves, associated 
to the transition. In Fig. [5][a) we plot these peaks for 
a polymer consisting of N = 100 monomers and a string 
with <jf = 3/2 for various values of the substrate adhesion 
strength ef. In Fig. E^b), we plot the transition tempera- 
tures corresponding to the peak positions depending on ef . 
We find an almost linear increase of the adsorption tem- 
perature with adsorption strength, a behavior which was 
in the same way observed in a recent study, where this 
transition has been studied for the same polymer model 
interacting with planar substrates Q. 

A more adequate description of the thermodynamic 
behavior of small finite-size systems provides the micro- 
canonical analysis 3, 13- III, based on the fact that all 
the relevant information about the system is encoded in 
its density of states g{E). In Fig. HJa), we plot the loga- 
rithm of this function, which is proportional to the micro- 
canonical entropy [S(E) = kslng^E)]. The adsorption 
transition is represented by the convex region in the micro- 
canonical entropy. The energetic width of this convex 
part is a measure for the latent heat which is nonzero for 
a first-order-like transition. The derivative of the micro- 
canonical entropy with respect to the energy yields the 
inverse microcanonical temperature (3(E) = dS(E)/dE. 
It is plotted in Fig. [6th) and exhibits a monotonic change 
in the transition region (therefore called "backbending ef- 
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Figure 5: (a) Adsorption peaks in the canonical heat capacities of 
the 100- mer for different values of tf at erf = 3/2 [cp. the peak 
temperatures with temperatures from the microcanonical analysis 
in Fig.[6lb)]. (b) "Phase diagram", i.e., transition line between ad- 
sorbed and desorbed phases. Constructed from peak positions in (a). 



feet" ) . The microcanonical transition temperature is ob- 
tained by a Maxwell construction in that region, indicated 
by horizontal lines in Fig. [6][b). We note, in particular, 
that for ef = 1 and 2, backbending does not occur. The 
inflection points in this energetic region indicate second- 
order-like transitions and belong to the O transition in 
bulk. This agrees with previous observations that the 
"strength" of the first-order signal decreases with decreas- 
ing substrate adsorption strength [34|- A more detailed 
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Figure 6: (a) Logarithm of the density of states (proportional to the microcanonical entropy) for the system with N = 100 monomers, 
crf = 3/2 and different string adsorption strengths ef. The convex intruder clearly emerges and becomes larger for increasing values of €f. 
(b) The derivatives of the functions in (a) with respect to E (proportional to the inverse microcanonical temperature). The lines mark the 
respective microcanonical adsorption temperatures obtained by the Maxwell construction in the backbending region. 



analysis of the structural transitions in the polymer-wire 
system by means of microcanonical thermodynamics will 
be presented in a forthcoming paper (lij . 
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